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We investigate the impact of an electric field on the structure of ultralong-range polar diatomic 
Rydberg molecules. Both the s-wave and p-wave interactions of the Rydberg electron and the neutral 
ground state atom are taken into account. In the presence of the electric field the angular degree 
of freedom between the electric field and the internuclear axis acquires vibrational character and 
we encounter two-dimensional oscillatory adiabatic potential energy surfaces with an antiparallel 
equilibrium configuration. The electric field allows to shift the corresponding potential wells in 
such a manner that the importance of the p-wave interaction can be controlled and the individual 
wells are energetically lowered at different rates. As a consequence the equilibrium configuration 
and corresponding energetically lowest well move to larger internuclear distances for increasing 
field strength. For strong fields the admixture of non-polar molecular Rydberg states leads to the 
possibility of exciting the large angular momentum polar states via two-photon processes from the 
ground state of the atom. The resulting properties of the electric dipole moment and the vibrational 
spectra are analyzed with varying field strength. 

PACS numbers: 31.50.-x, 33.20.Tp, 33.80.Rv 



I. INTRODUCTION 

Ultracold atomic and molecular few- and many-body 
systems offer a unique platform for a detailed under- 
standing and analysis of fundamental quantum proper- 
ties. The preparation and control of such systems in 
specific quantum states offer many opportunities for ex- 
ploring elementary quantum dynamical processes. Ex- 
perimentally, one can control the external motion of the 
atoms by designing and switching between almost arbi- 
trarily shaped traps [THS], and the strength of the in- 
teraction among the atoms can be tuned by magnetic 
or optical Feshbach resonances A striking new 

species are the weakly bound ultralong-range diatomic 
molecules composed of a ground state and a Rydberg 
atom whose existence has been predicted theoretically 
more than a decade ago and which have been dis- 
covered experimentally only recently [S]. The molec- 
ular Born-Oppenheimer potential energy curves, which 
are responsible for the atomic binding, show for these 
species a very unusual oscillatory behavior with many 
local minima. The latter can be understood intuitively 
and modeled correspondingly as the interaction of a neu- 
tral ground state atom with the Rydberg electron of the 
second atom. In a first approximation, the interaction 
between the two constituents is described by a s-wave 
scattering dominated Fermi-pseudopotential [H HQ] . The 
equilibrium distance for these molecular states is of the 
order of the size of the Rydberg atom and the vibrational 
binding energies are in the MHz to GHz regime for princi- 
pal quantum numbers n sa 30 — 40 depending on the type 
of states. More specifically, low-angular momentum non- 
polar states and large angular momentum polar states, 
so-called trilobites, have been predicted [1,. They pos- 
sess electric dipole moments in the range of IDebye (low- 
£)[TT] up to IkDebye (high-£)[7] in the polar case. The 



large electric dipole moment of the latter makes them 
accessible for electric field manipulation which opens the 
interesting possibility for the external control of molecu- 
lar degrees of freedom. Beyond the s-wave interactions, 
p-wave scattering has been shown [12, to lead to a class of 
shape-resonance-induced long-range molecular Rydberg 
states. 

The impact of magnetic fields on these ultralong-range 
molecules has been studied in |13| where it has in par- 
ticular been shown that the magnetic field provides an 
angular confinement turning a rotational degree of free- 
dom into a vibrational one and yields, with increas- 
ing strength, a monotonic lowering of the magnitude of 
the electric dipole moment. Polyatomic ultralong-range 
molecules formed of a Rydberg atom and several ground 
state perturber, such as coUinear triatomic species, can 
be constructed by taking the diatomic wave function as a 
basic unit and constructing the corresponding symmetry- 
adapted orbitals [TJ . Recently the formation of Rydberg 
trimers and excited dimers bound by internal quantum 
reflection [151 have been observed experimentally and an- 
alyzed in detail theoretically [T^. Moreover, it it has 
been shown how the electric field of a Rydberg electron 
can bind a polar molecule to form a giant ultralong-range 
stable triatomic molecule [l6l - fT8| which can consequently 
be controlled by applying external electric fields. Com- 
bining electric and magnetic fields in a crossed field con- 
figuration the existence and properties of so-called giant 
dipole ultralong-range molecules have been shown |19) . 
Opposite to the above-mentioned Rydberg molecules this 
species has no open radiative decay channels. 

In spite of the diversity of works focusing on ultralong- 
range molecules an original investigation of the impact of 
external electric fields specifically on the polar trilobite 
states is missing. Such an investigation is particularly 
desirable due to the strong sensitivity of these Rydberg 
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molecules with respect to the external field which pro- 
vides a handle on the control of their properties on a 
single molecule basis but also for their interactions in 
potential many-body systems. For these reasons we per- 
form in this work a study of the impact of an electric 
field on the structure and dynamics of high-^ ultralong- 
range molecules. We hereby proceed as follows. Section 
II provides a formulation of the problem presenting the 
working Hamiltonian and a discussion of the underly- 
ing interactions. Our analysis goes beyond the s-wave 
approximation and takes into account the next order p- 
wave term of the Fermi-pseudopotential. Section III and 

IV contain our methodology and a discussion of the ef- 
fects of the p-wecve contribution, respectively. In section 

V we analyze the evolution of the topology of the poten- 
tial energy surfaces (PES) with changing electric field. 
The resulting PES show a strongly oscillatory behavior 
with bound states in the MHz and GHz regime. With 
increasing field strength the diatomic molecular equilib- 
rium distance shifts substantially in a range of the order 
of thousand Bohr radii. We analyze the behavior of the 
corresponding electric dipole moment thereby achieving 
molecular states with a dipole moment up to several kDe- 
bye. Based on these properties and the s-wave admixture 
via the external electric field a preparation scheme for 
high-£ polar molecular electronic states via a two pho- 
ton excitation process is presented. Finally, we provide 
an analysis of the vibrational spectra which exhibit spac- 
ings of the order of several MHz. 

II. MOLECULAR HAMILTONIAN AND 
INTERACTIONS 

We consider a highly excited Rydberg atom interact- 
ing with a ground state neutral perturber atom (we will 
focus on the ^^Rb atom here) in a static and homoge- 
neous electric field. The Hamiltonian treating the Rb 
ionic core and the neutral perturber as point particles is 
given by (if not stated otherwise, atomic units will be 
used throughout) 

p2 

H= — +H,i + VnAr,n), (1) 

i/ei = ffo+Er, Ho = ^+Vi{r) (2) 

Where (M, P,R) denote the atomic Rb mass and the 
relative momentum and position of the neutral perturber 
with respect to the ionic core, (toc p, r) indicate the cor- 
responding quantities for the Rydberg electron. The elec- 
tronic Hamiltonian iJci consists of the field-free Hamil- 
tonian Hq and the usual Stark term of an electron in a 
static external E-field. Vi{r) is the angular momentum- 
dependent one-body pseudopotential felt by the valence 
electron when interacting with the ionic core. For low- 
lying angular momentum states the electron penetrates 
the finite ionic Rb'''-core which leads to a ^-dependence 
of the interaction potential Vi (r) due to polarization and 



scattering effects [20]. Throughout this work we choose 
the direction of the field to coincide with the z-axis of the 
coordinate system, i.e. E ~ Ee.^. Finally, the interatomic 
potential V^.c for the low-energy scattering between the 
Rydberg electron and the neutral perturber is described 
as a so-called Fermi-pseudopotential 

K,c(r, R) = 2T:A,[k{R)W " R) (3) 
-|-67rA3[fc(i?)]V(5(r-R)'^. (4) 

Here we consider the triplet (S* = 1) scattering of the elec- 
tron from the spin-^ ground state alkali atom. Suppres- 
sion of singlet scattering events can be achieved by an ap- 
propriate preparation of the initial atomic gas. In eq. Q 
As[k{R)] = -tan{5o{k))/k and Al[k] = - tan((5i(fc))/F 
denote the energy-dependent triplet s- and p-wave scat- 
tering lengths, respectively, which are evaluated from the 
corresponding phase shifts 5i{k), I = 0,1. The kinetic en- 
ergy i?kin — fc^/2 of the Rydberg electron at the collision 
point with the neutral perturber can be approximated 
according to fc^/2 = 1/R— l/2n^. The behavior of the 
energy-dependent phase shifts 5i as functions of the ki- 
netic energy -Ekin is shown in Fig. [l] 

III. METHODOLOGY 

In order to solve the eigenvalue problem associated 
with the Hamiltonian ([T]) we adopt an adiabatic ansatz 
for the electronic and heavy particle dynamics. We write 
the total wave function as \E'(r,R) — ■0(r; R)(/)(R) and 
obtain within the adiabatic approximation 

[Ho + Er + K.o(r, R)]V^.(r; R) = e,(R)7^,(r; R), (5) 

p2 

(— + e,(R))0,fc(R) = E,k<j)^k{'R), (6) 

where ipi describes the electronic molecular wave func- 
tion in the presence of the neutral perturber for a given 
relative position R and (j)ik determines the rovibrational 
state of the perturber. To calculate the potential en- 
ergy surface ei(R) we expand 7/;(r;R) in the eigenba- 

sis of Ho, i.e. ■0j(r;R) = J2nim ^nL(^)Xnini{r) with 

HoXnlrnir) £niX";m(l"), Xn/m(r) = Rnl{r)YimiO , (/)) . 

For / > /niin = 3 we neglect all quantum defects, i.e. Ho 
is identical to the hydrogen problem. Finally, we have to 
solve the following eigenvalue problem 

(£„/ - e(R))C„;„i -I- ^ Cn'l'rn'{E{nlm\z\n'l'm)5rr,m' 
n'V'in' 

+ {nlm\V^^^{Y,'R)\n'l'm')) = 0, (7) 

for which we use standard numerical techniques for the 
diagonalization of hermitian matrices. Throughout this 
work we mainly focus on the high-^ n — Zb manifold 
which provides the trilobite states in case of zero electric 
field To ensure convergence we vary the number of 
basis states finally achieving a relative accuracy of 10~^ 
for the energy. For the n = 35 trilobite manifold we used. 
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in addition to the degenerate n = 35, I > 3 manifold, a 
basis set that includes the 38s, 37d, 36p quantum defect 
split states due to their energetically closeness. This ba- 
sis set contains 1225 states in total. 
From eqs. ([5| and Q we already deduce some symme- 
try properties of the states tp, cf) and the energies e. If 
^r.R.E denotes the generalized parity operator that trans- 
forms (r, R, E) (-r, -R, -E) we have [i?, Pr,R,E] = 
[Vn,e(r, R), ^'r,R,E] = 0. This means that the states ^, tjj 
and (j) are parity (anti)symmetric and the PES fulfill 
e(R;E) = e(— R;— E). In addition, the PES possess an 
azimuthal symmetry, e.g. the vector defining the inter- 
nuclear axis can, without loss of generality, be chosen to 
lie in the x — z-plane. In the absence of any field, the 
PES depend exclusively on the internuclear distance R 
[7]. However, if a field is present, the PES are cylindrical 
symmetric, which means they also depend on the angle 
of inclination 9 between the field vector and the internu- 
clear axis, e. g. e(R) = e{R,6). 



IV. DISCUSSION OF THE P-WAVE 
CONTRIBUTION 
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FIG. 1. Energy dependent triplet phase shifts 5i and So for 
e~— *^Rb(5s) scattering. For — 24.7meV the phase shift 
Si — it/2, i.e. the energy dependent p-wave scattering length 
Ap{k) = — tan(5i(fc))/fc'^ possesses a resonance at this energy. 
This can be clearly seen in the inset, (color online) 



In several previous works the interaction between the 
Rydberg electron and the neutral perturber has been 
modeled to consist exclusively of a s-wave scattering po- 
tential [71[T3]. In this case, a single potential curve sphts 
away from the — degenerate high-/ manifold form- 
ing a strongly oscillating Born-Oppenheimer potential 
energy surface with a depth of around — 3.5 x 10^ GUz/n^ 
providing rovibrational states with a level spacing of ap- 
proximately lOOMHz and a permanent electric dipole mo- 
ment of IkDebye. However, it is an important fact that 
the p-wave scattering length Ap possesses a resonance at 
i?kin = 24.7meV. This corresponds to a radial distance of 
756ao for a Rydberg electron in a n = 35 state. The ef- 
fect of the p-wave contribution and its resonant behavior 
on adiabatic potential curves has been studied in detail 
in several works, see [12j[2T]. It has been shown that due 
to this interaction additional PES split away from the de- 
generate manifold leading to avoided crossings with the 
pure s-wave dominated potential curve and even ener- 
getically lower s, p and d-states. These avoided crossing 
cause a dramatic change in the topology of the PES lead- 
ing to novel molecular dynamics. 

According to [H] [32], for large radial distances the 
overall behavior of s and p-wave scattering dominated 
PES can be well described by the Borodin and Kazansky 
model [221 
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The energy curves are roughly determined by the corre- 
sponding phase shifts for the electron-neutral atom scat- 
tering, although the individual oscillations necessary for 
the existence of stable molecular states cannot be de- 
scribed by the simple equation ((sl). It can be shown that 



the positions of the crossing point of both PES Rc^^ and 
the position of the minimum of the s-wave dominated 
PES obey the ratio 



(n) 
c 



(m) 



(c) 



l/2n2 
l/2n2 



(9) 



with E^ 



,(c) 
kin 



7.1104nieV and E^'^J 



7.1102meV being 



the corresponding electronic kinetic energies (see Fig.[T]). 

Since 0.999973 < i?^"Vi?m^ < 1 it is in particular not 
possible to substantially shift the avoided crossings be- 
tween these two PES by increasing the principal quan- 
tum number n in the absence of the external field. This 
means that the p-wave contribution is crucial for the over- 
aU topology of the PES. 



V. POTENTIAL ENERGY SURFACES 

We focus here on the regime of field strengths E — 
— 650 — . The dissociation limits correspond to the 
atomic states Rb(5s)-|-Rb(n = 35,/ > 3) where / is used 
as a label in the presence of the electric field. In Fig. [2] we 
present the PES for the electrically dressed polar trilo- 
bite state for E = 150— and 300— as a function 9 and 
R. As mentioned in section HI in the absence of an elec- 
tric field the potential curves are independent of 9. For 
a finite field strength, this spherical symmetry is broken 
which is clearly seen in Fig. [2j For all field strengths 
the potential minimum is taken for the antiparallel field 
configuration 9 = tt. This is reasonable because the ex- 
ternal electric fields forces the electron density to align 
in its direction which leads to a higher density in the 
negative z-direction. The electric field therefore turns a 
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FIG. 2. Two-dimensional potential energy surfaces for the elec- 
trically dressed polar trilobite states for E — 150;^ (a) and 
300^ (b). We observe a potential minimum at 6 — tt. An 
increase of the electric field goes along with a stronger confine- 
ment of the angular motion and an increase of the diatomic 
equilibrium distance Rcq. Thus, the electric field stabilizes the 
s-wave dominated molecular states, (color online) 



rotational degree of freedom 6* to a vibrational one. As 
the field strength increases a stronger confinement of the 
angular motion is achieved and the corresponding equi- 
librium distance Rgq increases substantially. 
In Fig. ^) we show intersections through the PES for 
the 9th-15th excitation for a field strength of 300— for 
9 = TT. The inset in this figure shows the high-£ field- 
free potential curves. Due to the p-wave interaction 
a single potential curves splits away from the degener- 
ate manifold causing avoided crossings in the region of 
R = 1400ao — ISOOao. In addition, the inset also shows 
the potential curves provided by the Borodin-Kazansky 
approximation (see eq. ([s])). In the main figure the low- 
est potential curve is the one belonging to the 38s quan- 
tum defect split state. This state possesses a very weak 



oscillatory behavior in the MHz regime (not visible in 
Fig. (3^)). In general, this state is much less affected 
by the electric field compared to the PES arising from 
the zero field high-^ degenerate manifold. This is rea- 
sonable since the atomic 38s state does not possess a 
substantial electric dipole moment in the presence of the 
field. Therefore its potential curve hardly shifts with in- 
creasing electric field strength from its field-free value of 
— 20.36GIIz. However, the potential curves arising from 
the n = 35, / > 3 manifold show a strong dependence on 
the electric field. In general we obtain potential curves 
with a strongly oscillatory structure in the many hundred 
MHz to GHz regime. It is important to note that the p- 
wave interaction dominated PES (red curve in Fig. ([3^)) 
and the s-wave interaction dominated PES (blue curve in 
Fig. (|3^)) exhibit an avoided crossing which is crucial for 
the stability of the corresponding vibrational states. For 
zero field this avoided crossing happens to be compara- 
tively close to the global minimum of the corresponding 
'trilobite' PES (see inset in Fig. (|3^)) and might influence 
extended vibrational states. For increasing fleld strength 
however this avoided crossing is increasingly separated 
from the global equilibrium distance, as we shall discuss 
in more detail below. 

In Fig. ^p) we show intersections for 9 — tt through the 
PES of the electrically dressed polar states for different 
field strengths E = 150, 300, and 450^. We observe 
how the potential curve is globally shifted with increas- 
ing electric field strength. For E > 700^ (not shown in 
Figure) this trilobite PES experiences an avoided cross- 
ings with the potential curve belonging to the 38s state. 
In addition, we see that the overall topology of the PES 
do not change with varying i5-field. In particular, the 
number of minima and their positions remain approx- 
imately constant with increasing field strength. How- 
ever, the diatomic equilibrium distance i?eq (which is the 
global minimum in the range 1550ao < R < 2450ao) 
changes strongly as E varies. This means that the elec- 
tric field causes a 'spatial weight' to the PES. Low- 
lying vibrational molecular states in the well around the 
global minimum are shifted away from the region of the 
avoided crossings which is located in the interval 1500ao 
to ITOOoq. As a consequence the importance of the p- 
wave interaction is decreased significantly for higher field 
strength and the PES are determined mainly by the s- 
wave interaction. To be more specific we show the de- 
pendence of Rcq as a function of E in the inset in Fig. 
|4] We see a plateau-like structure with steps at the field 
strengths 100, 200, 385^ where the value of Rcq sharply 
changes. This structure simply reflects the depicted effect 
of the electric field on the PES, i.e. by varying the electric 
field one changes the energetically position of the differ- 
ent potential wells in the oscillating PES, which leads to 
abrupt changes of the global equilibrium position i?cq- 
Figs. ([3^,b) also demonstrate that with increasing field 
strength the avoided crossing between the p-wave inter- 
action dominated and the s-wave interaction dominated 
states remains (approximately) localized in coordinate 
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VI. ELECTRIC DIPOLE MOMENT 

In ref.[7] the authors reported on large electric dipole 
moments of ultralong-range polar Rydberg molecules of 
the order of kDebye. The zero-field permanent dipole 
moment for these species scales according to the semi- 
classical expression Dei = -Rcq — In Fig. |4]we show 
the absolute value of the electric dipole moment along 
the z-axis as a function of the field strength 

D,,{E) = I y"d3r7A*(r;Req,E)zVXr;Req,E)| 

C'n'l'mCnln, j drr^Rn'v {r)Rni {r) 

nn'll'ra 
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The integration over the angular degrees of freedom 
provides A/ = ±1 as a selection rule. We observe that 
with increasing electric field also Dei increases up to 
values of around 4kDebye. As for i?eq we see a sharp 
step structure, i.e. for field strengths at approximately 
100, 200, 385^ its values suddenly increase in steps 
of roughly 500 Debye. In Fig. |4] we also show a com- 
parison between the exact result calculated according 
to eq. ( 10 1 (blue data points) and the semi-classical 
approximation (green data points). For low electric 
fields the agreement is quite well, but differs with 
increasing field strength up to a deviation of around 
10%. The semiclassical result therefore certainly allows 
for a qualitative description of the behavior of D^x- 
For E > 570^ we find an unexpected decrease of 



FIG. 3. (a) Intersections through the two-dimensional PES for 
e = TT for the 9th-15th excitation for E = 300 — . For the 

m 

two lowest PES arising from the high-^ degenerate manifold a 
strongly oscillatory behavior is visible. The 38s split PES is only 
weakly oscillatory and hardly affected by the electric field. The 
inset shows the field-free trilobite and first p-wave PES. The re- 
gion of avoided crossings is clearly visible at _R ~ 1400 — ISOOao. 
In addition, we present the potential curves provided by the 
Borodin-Kazansky model given by eq. (|8|. (b) Same as in (a) 
but with varying E — '"^'^ '^"'^ /fsnV rj^\^„ 
librium distance R, 
avoided crossing at R 



150, 300 and 450^. The diatomic equi- 
cq is moving away from the region of the 



1500ao — 1700ao- (color online) 



space whereas the energetically low-lying potential wells 
with bound vibrational states and in particular the one 
belonging to the global equilibrium position are lowering 
in energy and are consequently well-separated from this 
avoided crossing. In conclusion, the electric field repre- 
sents an excellent tool to control the energetic positions 
and depths of the individual wells and to avoid destabi- 
lizing avoided crossing. 
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FIG. 4. The electric dipole moment as a function of the electric 
field E (blue points). For comparison we show a semi-classical 
prediction (green points). The inset shows the behavior of the 
equilibrium distance _Rcq with varying electric field strength, 
(color online) 

Del. This feature can be understood if one analyzes the 
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field-dependent spectrum of coefficients for tlie electronic 
eigenvector V(r;Req,E) = C'i(^)Xi(r)- In Fig. |6]we 
show the distribution |Cip for E = 300^ and 600^. 
For E ^ 300^ the spect rum is dominated by basis 
states from the n = 35, I > 3 manifold. Contributions 
stemming from the quantum defect split states (which 
are placed at the outermost right edge of the spectrum 
at i = 1217, ...,1225) are negligible. For E = 600^ the 
situation has changed in the sense that now the main 
contribution is provided by the 38s state. This can 
be understood by the fact that the considered PES is 
approaching the 38s PES with increasing field strength. 
The latter is however barely affected by the electric 
field. For E — 600— the PES involve avoided crossings 
which causes the high-£ dressed trilobite PES to acquire 
a major contribution from the 38s state. This finite 
admixture has two important consequences: 



Due to the A/ = ±1 selection rule the 38s state 



only acquires a contribution to the integral ( 10 ) 
via the 37p state. However, the coefficient of the 
latter state is negligibly small. This causes the 
decrease of Doi for large field strengths as seen in 
Fig.i 



• The finite 38s admixture provides us with the pos- 
sibility to prepare high-€ Rydberg molecules via a 
two-photon process. This goes beyond the three- 
photon preparation scheme suggested in [7] (^min = 
3 for the field-free case). For £; > 570^ the trilo- 
bite state acquires a major / = contribution 
which makes it accessible for a two-photon tran- 
sition scheme. The same mechanism has been re- 
ported recently in the analysis of ultralong-range 
polyatomic Rydberg molecules formed by a polar 
perturber [17]. Field-free high-£ molecular states 
can then in principle be accessed via an additional 
adiabatic switching of the electric field back to the 
zero value. 
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FIG. 6. Spectrum of coefficients of the electronic eigenvector 
tp{r; Req, E) at S = 300^ (a) and 600^ (b). For increasing field 
strength the eigenstates gain a finite admixture of the quantum 
defect split states. For 600^ we clearly see a major contribution 
provided by the 38s state, (color online) 



VII. ROVIBRATIONAL STATES 

Because of the azimuthal symmetry of the PES 
we introduce cylindric coordinates (p, Z, 0) for their 
parametrization e(R) = e{p,Z). For the rovibrational 
wavefunctions we choose the following ansatz 

(f>kumi'R-) = '^^"'"f^' ^•^ exp(zmy), m e Z, e Nq. 

Vp 

With this we can write the rovibrational Hamiltonian in 
eq. ^ as 



We solved the corresponding Schrodinger equation for 
different azimuthal quantum numbers m using a fourth 
order finite difference method. 

In Fig. [7] we provide the energies of the eleven lowest vi- 
brational (m = 0) states living in the trilobite PES for 
varying field strength. In order to obtain a normalized 
view of the spectrum the corresponding energy of the 
minimum of the PES has been subtracted. In general 
we observe a slight increase of the level spacing with in- 
creasing field strength. The increase is due to the an en- 
hanced angular confinement of the rovibrational motion 
for strong fields. For E = 100, 200 and 385— however we 
encounter a dip in the rovibrational level spacing. The 
latter corresponds to the case of crossover of the equilib- 
rium positions between neighboring wells and therefore 



an accompanying relocation of the corresponding rovi- 
brational wavef unctions. This leads to enhanced tun- 
neling probabilities between neighboring wells and there- 
fore an increased level density. In the inset of Fig. [7] we 
show the offset corrected potential curves for E — 300^ 
(blue curve) and E = 380^ (green curve) (9 — n). For 
E = 300^ the potential curve possesses a global mini- 
mum at i? = 1939ao and two local minima at i? = 1750ao 
and 2182ao with an offset of 200MHz. Bound states in 
the middle well with energies higher than 200MHz can 
tunnel into these wells, whereby their level spacing is re- 
duced. For increasing field strengths the right potential 
well is shifted downwards. This enhances the tunneling 
probabilities of states with energies less than 200MHz, 
which correspondingly leads to a denser spectrum. 
In Fig. (|8k,b) we present (scaled) probability densities 
\Fkvm{p,z)\'^ for m = for the vibrational ground state 
(v — 0) and the second excitation for E = 300 — . 
The equilibrium distance for the PES is located at 
Z = — 1939ao, p = 0. The to = ground state 
distribution is characterized by a deformed Gaussian 
profile that is localized a.t Z — — 1938ao and p — 
72ao . In Z, p-direction the density distribution pos- 
sesses an extension of approximately 50ao and lOOag, 
respectively. The density profile for the second excita- 
tion {v ~ 2) shows three separate Gaussian like den- 
sity peaks with increasing intensity located at {Z, p) — 
(-1939ao, 32ao), (-1934ao, 122ao) and (-1924ao, 234ao) 
with an extension of around (25ao, 30ao), (25ao,40ao) 
and (75ao,75ao) in the p-directions, respectively. 
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FIG. 7. Shown are the eleven lowest vibrational energies as a 
function of the field strength E. The dips around E = 100, 200 
and 385^ are caused by the change of potential wells deter- 
mining the diatomic equilibrium distance i?cq. In the inset we 
show the offset corrected potential curves for £ = 300-^ and 

^ m 

-B = 380-^ (9 = -k). For £ = 300-^ bound states in the mid- 
die well with energies larger than 200MHz can tunnel into the 
neighbored potential wells. This causes a reduction of their 
level spacings. For E = 380^ we nearly get a double potential 
well and states with energies less than 200MHz possess a higher 
tunneling probability. Correspondingly, this leads to a denser 
spectrum, (color online) 



VIII. CONCLUSIONS 

The recent spectacular experiments [HI [TTl [T5] prepar- 
ing, detecting and probing some of the important proper- 
ties of non-polar ultralong-range Rydberg molecules have 
opened the doorway towards a plethora of possibilities to 
create new exotic species where atoms, molecules or even 
clusters and mesoscopic quantum objects might be bound 
to electronic Rydberg systems. It is therefore of crucial 
importance to learn how the properties of these Rydberg 
molecules can be tuned finally leading to a control of the 
structure and potentially dynamics of these systems. The 
primary choice are here external fields in particular due 
to the susceptibility of the weakly bound Rydberg elec- 
trons. In the present work we have therefore explored 
the changes the polar high angular momentum trilobitc 
states experience if they are exposed to an electric field 
of varying strength. Taking into account s- and p-wave 
interactions it turns out that the electric field provides 
us with a unique know to control the topology of the adi- 
abatic potential energy surface. First of all, the angular 
degree of freedom between the electric field and internu- 
clear axis is converted from a rotational to a vibrational 
degree of freedom thereby rendering the field-free poten- 
tial energy curve into a two-dimensional potential energy 
surface. It turns out that the global equilibrium position 



is always the antiparallel configuration of these two axes. 
The sequence of potential wells with increasing radial 
coordinate, i.e. the oscillatory behavior of the potential 
is changed dramatically in the presence of the field. In 
particular we encounter an overall lowering of the en- 
ergy accompanied by a subsequently crossover of the en- 
ergetically order of the individual wells. Consequently, 
the equilibrium distance and the lowest vibrational states 
are systematically shifted to larger internuclear distances. 
The p-wave split state which, due to its resonant behav- 
ior, lowers dramatically in energy with decreasing inter- 
nuclear distance and therefore crosses the polar trilobite 
state close to its equilibrium distance in the zero-field 
case, can now with increasing field strength be systemat- 
ically shifted away from the energetically lowering equi- 
librium distance and corresponding well. In such a way 
the respective stability of the ground and many excited 
vibrational states of the polar trilobite state is guaran- 
teed. For strong fields the interaction of the latter state 
with non-polar (quantum defect split) states, which are 
very weakly polarized in the presence of the field, leads to 
a strong admixture of, in our specific case, s-wave char- 
acter to the polar high angular momentum states. As 
a consequence, a two-photon excitation process starting 
from the ground state of the two-atom system should be 
sufficient to efficiently excite these states and probe their 
character. The electric dipole moment, which is steadily 
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increasing with increasing electric field strength starting 
from zero-field, does, due to the above admixture, de- 
crease in the strong field regime. 

To obtain an even richer topology of the potential energy 
surfaces of the trilobite states the combination of static 
electric and magnetic field would be the next step. In this 
case both rotational degrees of freedom in field-free space 
will, in general, turn into vibrational modes rendering the 
field-free potential energy curve a three-dimensional po- 
tential energy surface. Depending on the configuration, 
such as parallel or crossed fields, the remaining symme- 
tries might even lead to controllable crossings or avoided 
crossings of the surfaces which will be the subject of a 
future investigation. 
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FIG. 8. Scaled probability densities \Fkvm{p, z)\'^ for rovibra- 
tional wavefunctions. Both wavefuntions belong to the trilobite 
PES for E = 300^ with an azimuthal quantum number m = 0. 
In (a) we observe a deformed Gaussian like density profile for 
the groundstate {1/ = 0) centered at Z — — 1938ao and p = 72ao. 
In Z/p-direction the density distribution has an extension of ap- 
proximately 50ao/100ao. In (b) we show the density profile for 
the second excited state (v = 2). This density profile provided 
three peaks at {Z,p) = (— 1939ao, 32ao), (— 1934ao, 122ao) and 
(— 1924ao, 234ao). (color online) 
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